(N 



Properties of ID classical and quantum Ising 

models: rigorous results 

o" 

^ ; William N. Yessen* 

^ ■ Department of Mathematics 

University of California, Irvine 

Irvine, CA 92697 

wyessen@math.uci.edu 



in 



Oh! 

■*— > 



> 



February 25, 2013 



Abstract 



In this paper we consider one-dimensional classical and quantum spin- 
1/2 quasiperiodic Ising chains, with two-valued nearest neighbor interaction 
modulated by a Fibonacci substitution sequence on two letters. In the quan- 
tum case, we investigate the energy spectrum of the Ising Hamiltonian, in 
presence of constant transverse magnetic field, by employing the techniques 
that were developed in [84]. In the classical case, we investigate and prove 
^sj ■ analyticity of the free energy function when the magnetic field, together 

£*\| with interaction strength couplings, is modulated by the same Fibonacci 

£S| ■ substitution (thus proving absence of phase transitions of any order at finite 

temperature). We also investigate the distribution of Lee- Yang zeros of the 
^— n | partition function in the complex magnetic field regime, and prove its Can- 

^vj . tor set structure (together with some additional qualitative properties), thus 

providing a rigorous justification for the observations in [3] and [5]. In both, 
quantum and classical models, we concentrate on the ferromagnetic class. 
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1 Introduction 

Since the discovery of quasicrystals [52,53,69,70], quasipcriodic models in math- 
ematical physics have formed an active area of research. The method of trace 
maps, originally introduced in [44,48,58] (see also [2,38,43,47,64,65,71,78] and 
references therein) , has provided a means for rigorous investigation into the phys- 
ical properties of one-dimensional quasiperiodic structures, leading, for example, 
to fundamental results in spectral theory of discrete Schrodinger operators and 
Ising models on one-dimensional quasiperiodic lattices (for Schrodinger operators: 
[11, 15, 17-21, 23, 48, 62, 72], for Ising models: [12, 13, 16, 24, 33, 75, 79, 85], and 
references therein). 

Quasiperiodicity and the associated trace map formalism is still an area of 
active investigation, mostly in connection with their applications in physics. In 
this paper we use these techniques to investigate crgodic families of classical and 
quantum quasiperiodic one-dimensional Ising models. 

1.1 One-dimensional quantum Ising chains 

One-dimensional quasiperiodic quantum Ising spin-1/2 chains have been investi- 
gated (analytically and numerically) over the past two decades [2, 12, 13, 16, 24, 
33,39-42,80,85]. Numerical and some analytic results suggest Cantor structure 
of the energy spectrum, with nonuniform local scaling (i.e. a multifractal) . The 
multifractal structure and fractal dimensions of the energy spectrum have only 
recently been rigorously investigated in [83], and only in some special cases. In 
this paper we extend the results of [83] to cover the general case. 

In the case of quantum Ising models, we investigate the energy spectrum of 
the Ising spin-1/2 chains with two- valued nearest neighbor couplings arranged 
in a quasiperiodic sequence, with uniform, transverse magnetic field. We shall 
concentrate on the quasipcriodic sequence generated by Fibonacci substitution on 
two symbols. 

By mapping the singlc-fcrmion Ising Hamiltonian canonically to a Jacobi ma- 
trix and employing the techniques that we developed in [84] , we rigorously demon- 



strate multifractal structure of the energy spectrum (thus extending the results 
of [83]-that is, answering [83, Conjecture 2.2]), and investigate its topological, 
measure-theoretic and fractal-dimensional properties. This provides analytic jus- 
tification for the results that were observed through numeric experiments (and, 
in some cases, rather soft analysis) by some authors - we note, in particular, 
[12,13,16,24,85]. 

In fact, we consider a family of Ising models where interaction strength cou- 
plings arc modulated by sequences from the hull of the Fibonacci substitution 
sequence, and by employing the results of [84], as mentioned above, we show that 
the spectrum is independent of the choice of a sequence in the hull. 

After investigating the energy spectrum of a single fermion, we pass to topo- 
logical description of the energy spectrum for many noninteracting fermions. In 
this case we show that for sufficiently many fermions, the spectrum is an interval 
whose length grows linearly as the number of particles is increased. 

1.2 One-dimensional classical Ising chains 

The classical Ising spin chain in one dimension with quasipcriodic interaction 
(nearest neighbor) and magnetic field have also been investigated both numeri- 
cally and analytically (see, for example, [2,3,5,27-29,33,75-77,79], and references 
therein). A particularly curious problem is analyticity of the free energy func- 
tion in case the interaction strength couplings as well as the external field are 
modulated by a quasiperiodic sequence (Fibonacci, say). Since the corresponding 
transfer matrices do not in general commute, this becomes a relatively nontrivial 
problem in comparison with the periodic case (we comment in more detail on this 
below, or see [2]). By employing the trace map as an analytic map on the three- 
dimensional complex manifold, wc show that the free energy function is in fact 
analytic in the thermodynamic limit. 

In fact, we consider a family of classical onc-dimcnsional Ising models where 
interaction strength couplings and magnetic field are modulated by a quasiperiodic 
sequence from the hull of the Fibonacci substitution sequence. It is well known 
that the partition function can be computed as the trace of the associated transfer 
matrices. These transfer matrices form a GL(2,R) cocycle over the shift map on 
the hull. Wc show that in the thermodynamic limit, the free energy is analytic, is 
independent of the choice of a sequence in the hull, and in fact coincides with the 
Lyapunov exponent of the associated GL(2, R) cocycle. 

Another problem is to describe topology of the Lee- Yang zeros of the partition 
function in the thermodynamic limit, when the former is viewed as a function of 
complex magnetic field. By the famous theorem of T. D. Lee and C. N. Yang 
[51], these zeros lie on S 1 C C when the couplings are ferromagnetic. Numerical 
experiments as well as some analysis in [3] suggest, among other things, that the 
distribution of zeros in the thermodynamic limit forms a Cantor set that exhibits 
gap structure similar to those that were investigated by J. Bellissard et. al. in [10] 
(see also [8,9,11] and [62]), when the sequence of interaction strength couplings is 
modulated by the Fibonacci substitution. This is also suggested by some results 
of J. C. A. Barata and P. S. Goldbaum in [5]. Both aforementioned investigations, 



however, fall short of proving rigorously the postulated gap structure or fractal 
properties of the Lee- Yang zeros. 

We prove that the Lee- Yang zeros in the thermodynamic limit do indeed form 
a Cantor set. In fact, we show that topological, measure-theoretic and fractal- 
dimensional properties of the set of Lee- Yang zeros on the circle are the same as 
those of the single-fcrmion energy spectra of quantum Ising quasicrystals, which 
are the same as those of the spectra of recently considered quasiperiodic Jacobi 
matrices. 

A natural question remains: are the Lee- Yang zeros in the thermodynamic limit 
independent of the choice of a sequence in the hull? At present we are unable to 
answer this definitely, but we conjecture a positive answer. 

2 The Fibonacci substitution sequence, symbolic 
dynamics on the hull and Lyapunov exponents 
of associated GL(k, R) cocycles 

Let us begin by constructing the Fibonacci substitution sequence, the hull of this 
sequence, and the associated symbolic dynamics (a left shift map on the hull). 
These constructions will play a crucial role in the rest of our investigation (to each 
sequence in the hull we shall associate an Ising model with interaction strength cou- 
plings and magnetic field modulated by the chosen sequence) . We then investigate 
Lyapunov exponents of GL(fc, R) cocycles over the shift map on the hull (this will 
play a role in our investigation of the free energy function in the thermodynamic 
limit of classical quasiperiodic Ising models). 

2.1 The Fibonacci substitution sequence and symbolic dy- 
namics on the hull 

Let A = {a, b}; A* denotes the set of finite words over A. The Fibonacci substi- 
tution S : A — > A* is defined by S : a H> ab, S : b i-> a. We formally extend the 
map S to A* and .4 N by 

S : aia 2 ■ ■ ■ a k H> S(cti)S(a 2 ) ■ ■ ■ S(a k ) 

and 

S : a\a 2 •••!-> S(u\)S(a.2) ■ ■ ■ ■ 

There exists a unique substitution sequence u S „4 N with the following properties 
[61]: 

Mi • • • up k = S 1 (a); S(u) = u; u\ ■■■ UF k+2 = U\ ■■■ UF k+1 U\ ■■■ UF k > (1) 

where {-Ffc} fceN is the sequence of Fibonacci numbers: Fq = F\ — 1; F k >2 = 
Fk-i + Fk-2- In fact, this sequence can be obtained constructively via the following 
recursive procedure. Starting with a, we apply the substitution s recursively: 

a i — > ab i — > aba t-^— > abaab i — y abaababa i — y ■ ■ ■ . 



From now on we reserve the notation u for this specific sequence. 

Let us remark that the sequence u can also be obtained as a sampling of an 
irrational circle rotation on two intervals as follows. Let (j> — (V5 + l)/2 denote 
the golden mean. Dchne the sequence u = {u n } nSN by 

u„ = ax[\-4,,i)(n(j) mod 1) + bx[i-<f,,i)c(n<f> mod 1), 

where » c denotes the set complement. 

Observe that for any Fibonacci number Ff~, k > 2, the finite word [u]i ■ ■ ■ [u]F k 
contains F^-i symbols a and Fk—2 symbols b. 

So far u is a sequence over N. Let us extend u arbitrarily to the left, and call 
the resulting sequence over Z, u. Endow A with discrete topology, and A with 
the corresponding product topology, and let T : A 1, <-^> denote the left shift. Define 
the hull of u, f2, as follows. 

ft = jw £ A z : lim T 7H (u) = uj; where {nj C N,m t oo} . 

Note also that there exists a special sequence cu s £ fl such that {[a) s ]„}„> 1 coin- 
cides with u. Indeed, 

ui s = lim T k (u). 

k— >oc 

From now on the notation lo s is reserved for this special sequence. 

It is easy to see that f2 C A % is compact and invariant under T. Moreover, 
it is known that T\q is strictly ergodic (that is, T\q is topologically minimal and 
uniquely ergodic) — for more details, see [61]. Let us reserve the notation \i for the 
unique ergodic measure for Ton!]. 

From now on whenever we make measure-theoretic statements about T on fi, 
such as, for example, existence of Lyapunov exponents almost surely, we implicitly 
assume that the statement is made with respect to the measure /i. We shall refer 
to the dynamical system T\q by (T, 0), and (T, f2,/x) when the measure fj, needs 
to be emphasized. 

2.2 Lyapunov exponents for GL(fc, R) cocycles over (T,Q) 

Define a measurable map G : Vt — >• GL(fc,M), such that for each u £ fl, G(uj) ^ 0. 
For n £ N, define g n : Q. — > R by 

g n {w) - log ||G(T"- 1 w)G(r' l - 2 a;) • • • G{Tlo)G{lo)\\ . (2) 

Observe that {g n }„ G N defines a subadditive process in the following sense. For 
any to £ Q, 

9n+m(w) < ffn(w) + g m {T n u). 

According to the quite general subadditive ergodic theorem of J. F. C. Kingman 
(see [46]), the limit 

lim -g n {u) 

n— >oo Ti 



exists /i-almost surely, is constant /i-almost surely, and belongs to [—00,00). In 
the literature this limit is often called the Lyapunov exponent of the GL(fe,K) 
cocycle over (T, 57,/i). By historical convention, the Lyapunov exponent is often 
denoted by 7; however, we reserve the notation 7 for other things, and denote the 
Lyapunov exponent by C: 

C := lim —g n (uj) /i-almost surely in uj. (3) 

n— >oo n 

In what follows, we shall be concerned with the situation where 

• G : 57^GL(2,M); 

• For every u> € 57, G(ui) has strictly positive entries; 

• G(w) depends only on the first term of the sequence ui, hence G is continuous. 

Then the following theorem by P. Walters (see [82]) guarantees existence and 
finitudc of the Lyapunov exponent not only almost surely, but everywhere. 

Theorem 2.1 (P. Walters, 1984). Let T : 57 — » 57 be a uniquely ergodic homeomor- 
phism of a compact metrizable space 57 and let G : 57 — > GL(fc,R) be continuous. 
If [G(w)]ij > for all u> € 57 and 1 < i,j < k, then C in (3) exists, is constant 
and finite for all u G 57. Moreover, if /x denotes the unique ergodic measure for 
(T,Q), then 

C= lim - / g n (uj)dfi(uj), (4) 

n— >oo n J 

where {g n } nt =m is as defined in (2). 

Let us mention quickly that according to A. Hof 's arguments [37, Proposition 
5.2], for substitution dynamical systems, such as (T, 57), generated by primitive and 
invertable substitutions, such as the Fibonacci substitution — see [61] for definitions 
and details — C exists and is constant for all w S 57 even when G(oj) contains zero 
entries; as long as there exists k > 1 such that G(ui) depends only on the finitely 
many terms {w_fc, . . . , ljq, . . . a;*;} of the sequence to. 

3 Description of the models and main results 

We begin with discussion of the quantum Ising model on a one-dimensional integer 
lattice, with constant transverse magnetic field and interaction couplings modu- 
lated by a quasiperiodic sequence chosen freely from the hull 57 of the Fibonacci 
substitution sequence u, as constructed in the previous section. 

3.1 One-dimensional quantum Fibonacci Ising model 

Let p : A —$■ R, with p(a),p(b) > and, unless stated otherwise, p(a) ^ p(b). For 
w G 57, set J Uyn = p([T n u>]i), and define a quantum Ising Hamiltonian on a lattice 



of size N by 

JV N 

W^-^V^VW-ft^aW, (5) 

n— 1 n— 1 

where /i > is the external magnetic field in the direction transversal to the lattice. 
The mati 
given by 



The matrices ayf are spin-f /2 operators in the x and z direction, respectively, 



a (x)M =i(g)...<g,ig) o-(»).W (g) I® ... ®I , 

N—n terms nth position N-n+1 terms 

where I is the 2x2 identity matrix. Here qW'W arc the Pauli matrices given by 

...» = (; j) md .« . (j _« t 

that we view as operators on C 2 . Thus H u is a linear operator on <g>-product of 
N copies of C 2 . We deliberately do not insist on a specific choice of boundary 
conditions, since taking the thermodynamic limit (described below) renders any 
particular choice of boundary conditions redundant. 

Scaling by 1/h, by abuse of notation we can write J u . n = J u ,n/h, so that the 
Hamiltonian in (5) can be written as 

N 

u » = - E (-w^Bi + <*$) • ( 6 ) 

The Hamiltonian in (6) can then be extended periodically to the infinite periodic 
lattice with unit cell \io\\ ■ ■ ■ [lo]n- Let us denote the resulting Hamiltonian by 
rLu,N- 

The energy spectrum of liu^N can be computed as the set of those values A 
that satisfy the equations 

(A-B)4> = \^ (7) 

{A + B)tJj = \<f>, 

where A and B are infinite tridiagonal periodic matrices given by 

(all other entries being zero) for 1 < i < N, and Ai+jf,j+N = Ai,ji ^i+TV.j+w = 
Bij. and the corresponding wave function (4>,ip) = {(4> n ,ifi n ) G C 2 } pN does not 
diverge exponentially (for more details see, for example, [24,54]). 
Observe that (7) is equivalent to 

(A-B)(A + B)i> = \ 2 4K 



This equation can be written as 

Ju,n-l1pn-l + (1 + Jl.Mn + J u ,n+l^n+l = A 2 t/>„, 

where J Wj j = J^j for 1 < i < iV, and J u ,i+N = Jw,*- Define an operator JCi by 

[ l 7i Af) V']« = Jw.n-lV'n-l + (1 + Jl,n)^TL + Juj, n+l^n+1- 

The operator JTLj is an infinite periodic Jacobi operator of period N. and as 
N — > oo, J7Lj converges in strong sense to the operator J u defined by 

[Ju1p]n = Ju,n-l^u,n-l + (1 + Jri)4'n + Jcuji+l4'n+l (8) 

(for more details see [84]). Consequently, the spectrum oiH^.N can be investigated 
as the spectrum of Ju , which, by Floquet theory, consists of N compact intervals 
(see, for example, [74]). Let us denote the spectrum of H^.n by q u ,n- We are 
interested in the structure of the spectrum in the thermodynamic limit (that is, 
as N — > oo). 

Before stating one of our main results, let us quickly recall definition of the 
Hausdorff metric, as a metric on 2 R and the definition of local Hausdorff dimension. 

For any A, B C K, define the Hausdorff metric distn(^4, B) by 

distn(^, B) = max < sup inf {\a — b\} , sup inf {\a — b\} > . 

laeAbeB beB aeA J 

Based on the connection mentioned above between <7 Wi at and the spectrum of 
3u , one may postulate a connection between the spectrum of J^ in (8) and the 
spectrum of % u ,n in the thermodynamic limit. This indeed follows, since the Ising 
model is equivalent, via a unitary transformation, to the tight binding model in 
(8) (see, for example, [50,67,68] for an attempt to approach Ising models via these 
tight binding models — mostly numerical). 

It was shown in [84] that the spectrum of J u is independent of u> E il. The same 
independence result follows for Huj.n in the thermodynamic limit. It is therefore 
enough to study T-L Us .n, where u s £ fl is as defined in Section 2.1. 

Denote by er^ the spectrum of - H tlJsi F fc , where F^ is the fcth Fibonacci number. 
It was shown in [84] that as k — > oo, the spectra of J^i, approach in Hausdorff 
metric the spectrum of J uls . Again, via the aforementioned connection between 
the Ising Hamiltonian and Jacobi matrices, one would expect a similar result for 
the Ising Hamiltonian. This is indeed the case. 

To avoid a rather opaque technical thicket, we shall proceed directly by explic- 
itly constructing a compact set and showing that the sequence {cfc} approaches this 
set in the Hausdorff metric. We then proceed to describe topological, measure- 
theoretic and fractal-dimensional properties of this limit set, which we call the 
energy spectrum of T-L Ub ,N in the thermodynamic limit. 

Before we continue, let us introduce the following notation. We shall denote 
the Hausdorff dimension of a set A C M by dim^A). By dixn^ c {A^a) we denote 
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the local Hausdorff dimension of A at a point a: 

diniH C ( J 4, a) — lim dim H (^4 n (a - e, a + e)). 

The box-counting dimension of A is denoted by dime (.A). 

Remark 3.1. In what follows, by a Cantor set we mean a (nonempty) compact set 
with no isolated points whose complement is dense. 

Since a certain multifractal with specific properties will appear more than once, 
it is convenient to give this set a name: 

Definition 3.2. We call a Cantor set Seta um-set (for uniform multifractal 
set) provided that S satisfies the following properties. 

(a) The set S has zero Lebesgue measure; 

(b) The function d that maps s £ S to dim I ^ c (S', s) is continuous and for all e > 
and for s G S, d restricted to (s — e, s + e) PI S is nonconstant (in particular, S 
is "uniformly" — i.e. localized at every s £ S — a multifractal); 

(c) The Hausdorff dimension of S is strictly between zero and one. 

If (c) fails, we call S full um-set. If S satisfies (a)-(c) and, in addition, 

(d) The box counting dimension of S exists and coincides with its Hausdorff di- 
mension, 

then we call S tamed um-set. If S satisfies (a)-(d) except (c) , we call S a full 
tamed um-set. 

Let us recall that, roughly speaking, a dynamically defined Cantor set C is one 
which is a blowup of an arbitrarily small neighborhood of any a € C under a C 1+e 
map defined in a neighborhood of C (here C 1+e means continuously diffcrcntiablc 
with Holder continuous derivative whose Holder exponent is e). For the technical 
definition and properties of dynamically defined Cantor sets, see [59, Chapter 4]. 
Let us only mention that dynamically defined Cantor sets are rather rigid. For 
example, box-counting dimension exists and coincides with the Hausdorff dimen- 
sion. Moreover, local box-counting and Hausdorff dimensions do not depend on 
the point of localization and coincide with the global dimensions. 

It is convenient to give a name to a Cantor set which is, in a sense, arbitrarily 
close to being dynamically defined: 

Definition 3.3. We call a Cantor set S almost dynamically defined Cantor set 
provided that there exists a finite subset T of S (possibly empty) such that outside 
an arbitrarily small neighborhood of 7", the rest of S is a dynamically defined 
Cantor set. 

We are now ready to state our main results for the quantum Ising quasicrystal. 



Theorem 3.4. There exists a compact (nonempty) set Boo in K such that the 
sequence {<7fc} feGN converges to Boo in the Hausdorff metric. The set Boo implicitly 
depends on the choice of p(a) andp(b), and has the following properties. 

i. .Boo is a um-set, but not full um-set; 

ii. With p(a) fixed, for all choices of p(b) sufficiently close to p{a), but not equal 
to p{a), Brx, is a- tamed um-set, but not full tamed um-set; 

Hi. dimH(-Boo) depends continuously on (p(a),p(b)). 

Remark 3.5. Thus Theorem 3.4 is an extension of [83, Theorem 2.1]. Also, we 
believe the restriction on parameter p{b) in statement ii can be dropped (see the 
concluding remarks in [84] for more details). 

Next we discuss topological structure of the energy spectrum for many nonin- 
teracting fermions. This can be thought of as TV identical Ising models on separate 
ID lattices with no interaction between them. For N > 1 noninteracting fermions, 
the corresponding energy spectrum is the TV-fold arithmetic sum of single-fermion 
spectra, minus the diagonal; that is, denoting by B^ the iV-fermion spectrum, 
we have 



Boa = i a i + • • • + a N '■ a-i € Boo and for any i ^ j, ai ^ a,}, 

where • denotes topological closure. On the other hand, since, as claimed in 
Theorem 3.4, Boo is a Cantor set, it contains no isolated points. Therefore, 

B H = Wi H h ajv : Oj € -Boo} ■ 

Regarding topological structure of B^ , we have the following. 

Theorem 3.6. For every choice of positive and distinct p{a), p(b), there exists 
No € N such that for all N > Nq, B^ is the interval [A, B], where 

A = min < ^ a^, a^ <E Boo > B = max < ^ b Ll b % € Boo \ ■ 

Consequently, for N > N , 

Blo = l-Nb,Nb], (9) 

where b is the upper bound of Boo! hence if N = oo, then B^ = K. Moreover, for 
any p(a) > and all p(b) > sufficiently close to p(a) and not equal to p(a), we 
can take Nq = 2. 

Remark 3.7. Topological structure of B^ for N < Nq is an interesting question. 
At present we do not have any strong qualitative results in this direction. We 
believe that when p(a) and p(b) are sufficiently far, then for some N > 2, B^ 
is a Cantorval — a compact set with dense interior, any connected component of 
the complement of which has its endpoints approximated by the endpoints of 
other connected components of the complement. In other words, it is in a sense 
a "Cantor set with dense interior." Sets of this nature, and their appearence in 
dynamical systems and number theory, have been under active investigation. 
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3.2 One-dimensional classical Fibonacci Ising model 

We first consider the classical ID Ising spin chain in real magnetic field and real 
temperature regime. For the sake of brevity, we refer to it simply as the real 
regime, and in the complex magnetic field case - complex regime. 

3.2.1 Real regime 

Let us consider a linear chain of N spins ui € {±1} , 1 < i < N . The energy of a 
configuration a = (<ti, a%, . . . , <7at) is given by 



N 

E 



H(a) = - V ( JiOiO-i+i + HiCJi) , (10) 



with circular boundary conditions (i.e. cttv+i = o~i)- As mentioned in the in- 
troduction, we concentrate our attention on the ferromagnetic case (that is, the 
interaction strength couplings Ji are positive). We take Hi, the external magnetic 
field, nonncgative. 

The partition function of this system is given by 



Z(-)=£ex P (_J-W), (ID 



k B r 

where the sum is taken over all possible configurations of N spins. Here r denotes 
the temperature and ks - Boltzmann's constant. The corresponding free energy 
is given by 

p(">= 1-logZW (12) 

Nkbt 

expressed as a function of r. 

Let p : A — > M, p(a),p(b) > and p(a) ^ p(b) unless stated otherwise. Also, 
introduce q : A — > R with q(a),q(b) > (we do not force, in general, that q{a) ^ 
q(b)). For any u) € O, For n G N, set 

J u , n =p([T n u]i) and H u , n = q([T n u] 1 ). 

With these sequences of interaction strength couplings, { J w , n }, and magnetic field, 
{-H^.n}, we obtain a classical (i.e. nonquantum) quasipcriodic Ising model in (10), 
which depends on the choice of u> € Q. To emphasize this dependence, let us 
write Zi, and Fi, for the corresponding partition function and the free energy, 
respectively, on a lattice of size N. 

Let us remark that the partition function Z^ N ' is obviously strictly positive, 
so that F^ is well-defined. Regularity properties of Ft"' in the thermodynamic 
limit (that is, as N — > oo), which we denote by F (or by F u when dependence on 
lo G Q needs to be emphasized), presents a curious problem in connection with 
critical phenomena (i.e. phase transitions). Since Z' w ' does not admit zeros for 
finite N and < r < oo, the free energy does not experience critical behavior in 
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finite volume and finite temperature (this is of course expected in dimension one) . 
In fact, for finite TV and < r < oo, F^ is easily seen to be analytic. Does 
this property survive the thermodynamic limit? It is expected that the answer is 
yes. In some cases F can be computed analytically; in the present case, however, 
properties of F^ (assuming F^ is well-defined) were hitherto not known (see, for 
example, remarks in [2]). This brings us to 

Theorem 3.8. The sequence of functions Ft, converges, as TV — > oo, uniformly 
in t to an analytic function F of t £ (0, oo) which is independent of u> and is 
strictly negative. 

Remark 3.9. Our proof of Theorem 3.8 also provides a recursive algorithm for 
numerical approximation of F. 

3.2.2 Complex regime 

In this case we take q : A — > C constant — denote this constant by H — and < 
t < oo. For a choice of a; £ f2, the partition function Zjj can then be considered 
as a function of three variables: 

•--»($)■ >"*($)■ '-«»(£) < 13 » 

(that is, a, f3 and 77 are functions of r: a, ft : r 1— > R, 77 : r i-> C). From now on we 
shall concentrate on fixed r £ (0,oo) and variable H £ C (that is, a,/3 are fixed, 
while the fugacity variable r\ is varied over C). 

While Till does not admit zeros in the real regime, it does (even in the fi- 
nite volume TV and finite r £ (0,oo)) in the complex regime, as a function of 77 
(indeed, ZL ' (rj) is a polynomial in rj). These so-called Lee- Yang zeros [51] are 
guaranteed to lie on the unit circle S 1 C C if both couplings are ferromagnetic 
(that is, p(a),p(b) > 0, which is the case considered here); similarly, in the purely 
antifcrromagnetic case (that is, p(a),p(b) < 0), the zeros lie on the negative real 
axis; the mixed case is much more subtle and is not treated here (for further details 
and a more general exposition see, for example, [66]). We are interested here in 
distribution of zeros of ZL in the thermodynamic limit, N — >• 00. 

Let Z^^n denote the set of zeros of Z^ in the variable rj. Some numerical 
experiments and soft analysis in [3] suggest that in the limit TV — > 00, {-Z Ws! jv} 
accumulates in a Cantor set on S , where uj s £ Q, is as defined in Section 2.1. 
Also, a gap structure similar to that studied by J. Bellissard et. al. in [10] 
was observed. We rigorously justify these observations, and add further detail to 
qualitative description of the zeros in the thermodynamic limit in the following 

Theorem 3.10. There exists a compact, nonempty set Z C S 1 C C such that 
Zuj e n 5- Z in Hausdorff metric. Moreover, the set Z has the following prop- 

N— s-oo 

erties. 
i. Z is a um-set together with (possibly) finitely many isolated points; 
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ii. {±1} n Z = (as expected, of course); 

Hi. With p{a) fixed, for all p{b) sufficiently close to p{a), Z is a tamed um-set; 

Remark 3.11. We believe that in statement i, there do not exist the mentioned 
isolated points, though at present we are unable to prove it rigorously; however, iii 
states that no isolated points are present provided that p(a)/p(b) ~ 1 (assuming 
P (b)^0). 

Since in Theorem 3.10 statements are made only for the special choice u> s E il, 
a natural question still remains: how is the conclusion of Theorem 3.10 affected 
as one considers general uj £ f2? At present we do not have a general answer, but 
we conjecture that the conclusion of the theorem is not affected by any specific 
choice of uj £ f2, and, moreover, Z in Theorem 3.10 is independent of the choice 
of u € CI. 

4 Proofs of main results 

For proofs of main results we rely heavily on dynamics of the so-called Fibonacci 
trace map (to be defined momentarily), and the geometry of some dynamical 
invariants. We recall here only those properties of the aforementioned objects 
that are necessary to prove our results, and only briefly. For a survey on trace 
maps (including the Fibonacci trace map) associated to substitution sequences, 
the reader may consult, for example, [4,64,65]. For the necessary notions from 
hyperbolic and partially hyperbolic dynamics (including notation and terminol- 
ogy), see a brief outline in Appendix B, and references therein for more details. 
For model-independent results (some of which appear below), see [22, Section 2]. 

4.1 Dynamics of the Fibonacci trace map: a brief outline 

Define the so-called Fibonacci trace map on C 3 by 

f(x,y,z) = (2xy-z,x,y). (14) 

Define also the so-called Fricke- Vogt character (often also called Fricke- Vogt in- 
variant) [25,26,81] I : C 3 ->• C by 

I(x,y, z) = x 2 + y 2 + z 2 -2xyz- 1. (15) 

It turns out that / is invariant under the action of /, that is, 

I(x,y,z) = Io f(x,y,z) 

(in connection with Schrodinger operators, see the pioneering work in [44,45,48, 
49,58]). Since / is invertible, it follows that the sets 

SveC = { (x, y, z) e C 3 : I(x, y, z) = V) 

are also invariant under /, and in fact f(Sy) = Sy- The sets Sy are complex- 
analytic surfaces. 
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In what follows, we sometimes consider the real part of Sy separately (that 
is, the subset Sy fll 3 of SV). Moreover, as far as the real part, we only need to 
consider S v n R 3 for V G R and V > 0. In this case, if V > 0, then S v n R 3 
is a smooth, connected, unbounded two-dimensional submanifold of R 3 without 
boundary. It is topologically the four-punctured 2-sphere. If V = 0, then Sv H R 3 
has four conic singularities. Minus these singularities, it is a smooth manifold with 
five connected components: one bounded, diffcomorphic to the four-punctured 2- 
sphere, and four unbounded, diffeomorphic to the once punctured 2-disc. 

To distinguish between the "full complex surface" Sy and Sy nR 3 , we denote 
the former by S v and the latter by Sy. We call Sq 1 the Cayley cubic. 

It isn't difficult to see that S v is also invariant under /. In what follows, we are 
concerned mostly with those points in Sy, that have bounded forward semi-orbit 
under the action of /; that is, p G S v such that 

OJ(p) := {f n ( P ) : n G N} 

is bounded (we define the backward semi-orbit ofp, OJ (p), similarly by considering 
iterates of / _1 in place of iterates of /). Following convention from the theory of 
holomorphic dynamical systems, let us define the following sets: 

Ky- := Ipe S v : Oj is bounded} and Ky := ip G Sy : OJ is bounded} . 

Define also 

Qv—K+nKy. 

It turns out that Kp G S v (see [14]). It follows that Oy C Sy. Also, if a point 
has unbounded forward semi-orbit, then the orbit escapes to infinity, i.e. the orbit 
does not contain a bounded subsequence (see, for example, [63]). 

Regarding the topological structure of fV and dynamics of / restricted to ily 
(clearly VLy is /-invariant) we have 

Theorem 4.1 (M. Casdagli [15], D. Damanik and A. Gorodetski [20], S. Cantat 
[14]). l For all V > 0, fly C Sy is a compact Cantor set of zero Lebesgue mea- 
sure. It is precisely the nonwandering set for f restricted to Sy and f on ily is 
hyperbolic. Moreover, fly is locally maximal and f\a v is topologically transitive. 

From [83] we have the following 

Proposition 4.2 ([83, Combination of Corollary 3.8 and Proposition 3.9]). For 
any V 2 >V\> 0, Uve[Vi,v 2 ] ^ v is V^rtially hyperbolic for f . Thus, Uve(Vi,V 2 ) K v 
is a family of pairwise disjoint smooth open two-dimensional injectively immersed 
submanifolds of R 3 without boundary, called the center-stable manifolds. The 
center-stable manifolds intersect the surfaces Sy, V > 0, transver sally. 



M. Casdagli proved this theorem for V > sufficiently large, D. Damanik & A. Gorodetski 
extended Casdagli's result to all V > sufficiently small, and S. Cantat proved the result for all 
V > 0. M. Casdagli and D. Damanik & A. Gorodetski used Alekseev's cone technique to prove 
hyperbolicity, while S. Cantat used other techniques from holomorphic dynamics. 
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The following result, while easily obtainable, is convenient to state in the form 
of a proposition for later reference. For a proof, see, for example, [84, Lemma 3.4]. 

Proposition 4.3. For oo > Vi > Vo > 0, the set Uy e ry vil ^v * s a c ^ ose d subset 
o/R 3 , and consequently also o/C 3 . 

We shall also need regularity of escape rate of those points whose forward 
semi-orbit is unbounded. 

Proposition 4.4. Let D be a domain in C 3 such that no point of D has bounded 
forward semi-orbit under f . Then the map £ : D — > R defined by 

log|l/"(z)|| f . 

£ : z ^ lim — , (16) 

n— >oo m 11 

where <f> = (vo + l)/2 is the golden mean, is well-defined and analytic. Moreover, 
£{z) = £(z), where 

^„ llra fe^i!M, (17) 

where it denotes projection onto the first coordinate. We assume that irof n (z) ^ 0, 
which is guaranteed for sufficiently large n by hypothesis. 

Proof of Proposition 4-4- Existence of the limit in (16) follows from [14, Corollary 
3.3]. That the sequence log ||/ n (^)|| /(f> n is uniformly Cauchy in z, and hence 
converges uniformly in z, follows from arguments similar to those given in the 
proof of [7, Lemma 3.2]. Since for every n > 1, the map z i-> log ||/ n (z)|| /0" i s 
pluriharmonic on D, it follows that its uniform limit, £, is also pluriharmonic, and 
hence analytic, on D. 

The second assertion (that is, (17)) follows from the fact that escape rate of 
orbits under / is the same as escape rate in each coordinate (for details, see, for 
example, [14,63]). 

□ 

Next we state one of the main results of [84] in a slightly modified, more general 
geometric form. This result will play a central role in the proof of Theorem 3.10 
below. Although proofs in [84] can be easily extended to cover this more general 
case, we invite the reader to see a discussion on model-independent results in 
[22, Section 2] for technical details. 

Theorem 4.5 ([84, Theorem 2.3]). Suppose 7 : R — > {J v>0 Sy is a compact 
regular analytic curve that intersects the center-stable manifolds. Then either 7 
lies entirely on a center-stable manifold, or it intersects the center-stable manifolds 
transversally in all but, possibly, finitely many points. Denote the set of tangential 
intersections byT. LetT denote the set ^C\{center- stable manifolds}. IfT contains 
points outside of T other than just the endpoints of 7, then the following hold. 

i. T is a um-set together with (possibly) finitely many isolated points if and only 
if 7 does not lie entirely in some Sy, for some V > 0; otherwise T is an 
almost dynamically defined Cantor set; 
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ii. If T = 0, then V is a tamed um-set if T does not lie entirely on some Sy, 
V > 0, and a dynamically defined Cantor set otherwise; 

Hi. r \ T is a full um-set if and only if there exists a point inT \T lying on Sq . 

iv. Suppose 7 depends analytically on a parameter k. Write T K for the corre- 
sponding set of intersections with the center-stable manifolds. If for some kq, 
r ' Ko is a um-set together with (possibly) finitely many isolated points, then 
for all k sufficiently close to kq, T k is also a um-set together with (possibly) 
finitely many isolated points. Moreover, ifT Ko does not contain isolated points, 
then for all k sufficiently close to kq, T Ka does not contain isolated points and 
K i— > diniH(r K ) is continuous. However, the properties of being a full um-set 
and almost dynamically defined Cantor set do not survive arbitrarily small 
K-perturbations. 

We are now equipped with all the tools needed to prove the main results. 

4.2 Proof of Theorem 3.4 

For any choice p(a),p(b) > (as usual, p(a) ^ p(b)), define the line 7 in R 3 
parameterized by E € M. by 

, E 2 -(l+ P (br) E 2 -(l+p(a) 2 ) p{af+p{b)< 
j{E) - 



Mb) Ma) ' 2p(a)p(b) 

Define also 

Boo = {e E E : (D+(~f(E)) is bounded | . 

It was shown in [83] that B m is compact (and, of course, nonempty) and that 
(Jk > Boo in Hausdorff metric. Thus the first claim of Theorem 4.2 follows. We 

k— >oo 

now investigate properties of Boo ■ We start by noting that -Boo contains infinitely 
many points none of which are isolated (this can be seen from relation of Boo 
with the spectrum of the tight-binding Hamiltonian in (8), which was investigated 
rigorously in [84] ) . 

Remark 4.6. Let us remark that properties of Boo were investigated in [83] in the 
special case that for fixed p(a) > 0, p(b)/p(a) is sufficiently close to one. 

It is of course enough to consider 

~ (j g )= (E-(l+p{bf) E-{l+p{af) p(af+p(bf\ 



Mb) Ma) ' 2p(a)p(b) J 

in place of 7, with E <G [0,oo). Evaluating the Fricke-Vogt invariant from (15) 
along 7, we get 

J (7$)) = f (;-rV, (18) 
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where r =p{a)/p{b). Clearly I(j(E)) <E R and I{^{E)) > 0. Moreover, we have 

nm) = U 1 --X>o. 



dE 4 \r 

Since the image of 7 coincides with that of 7, it follows that 7 is transversal to 
every S v , V > 0, and from (18) it follows that 7 intersects each S v , V > in only 
one point. 

The points of intersection of 7 with the center-stable manifolds lie on a compact 

line segment along 7, since for all E such that E is sufficiently large, j(E) satisfies 

the hypothesis of [83, Proposition 3.1-(2)], and hence escapes. We are therefore 
working in the setting of Theorem 4.5. 

Observe that 7 intersects Sq at E = 0, in the point 

With (xo,ho,zq) = f 2 (po), a simple (albeit long) computation shows |xo|,|j/o| 
and \z \ are greater than one, and \xoyo\ > \z \. It follows that Ot(/ 2 (p )) 
is unbounded (see [83, Proposition 3.1-(2)]). Thus Theorem 3.4-i follows from 
Theorem 4.5-iii. 

In [83] it was shown that with p(a) > fixed, for all p(b) > sufficiently close 
to p(a), 7 intersects the center-stable manifolds transvcrsally. Hence Theorem 
3.4-ii follows from Theorem 4.5-ii. 

Again, Theorem 3.4-iii follows from Theorem 4.5, since 7 depends jointly ana- 
lytically on p(a) and p(b). 

This completes the proof of Theorem 3.4. 

4.3 Proof of Theorem 3.6 

First we need to introduce the notion of Newhouse thickness (see [57]). 

Let C C R be a Cantor set and denote by 1 its convex hull. Any connected 
component of X\C is called a gap of C. A presentation of C is given by an ordering 
U = {U n } n >i of the gaps of C. If u € C is a boundary point of a gap of U of C, we 
denote by K the connected component of X \ {JJy U U% U • • • U U n ) (with n chosen 
so that U„ = U) that contains u, and write 



r(C,W,«) = & 
With this notation, the thickness t(C) of C is given by 



t(C) = supinf t(C,U, u) 
u u 
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Next we define what S. Astels in [1] calls normalized thickness. For a Cantor 
set C, define the normalized thickness of C by 

m- T(C) 



1 + t(C) 

(Astels originally used the letter 7 for normalized thickness, which we use here for 
the line of initial conditions to preserve notation from [83,84] for easy reference). 
A specialized version of the more general [1, Theorem 2.4] gives the following 
remarkable result. 

Theorem 4.7. Given a Cantor set C, if mv(C) > 1, m <G N, then the m-fold 
arithmetic sum of C with itself is equal to the m-fold sum of I with itself where X 
is the convex hull of C. 

In order to apply Theorem 4.7 in our case, it is enough to show that for any 
positive and distinct p(a), p(b), we have t{B oc ) > 0. 
For a fixed a € -Boo define the local thickness at a by 

lim t(Boo n (a — e, a + e)) 

(that the limit exists follows easily from the fact that thickness is decreasing with 
respect to set inclusion). From [84, Lemma 3.6] it follows that for any a £ Boo, the 
local thickness at a coincides with thickness of some dynamically defined Cantor 
set. On the other hand, thickness of any dynamically defined Cantor set is strictly 
positive (see [59, Chapter 4]). Hence t(B 00 ) > 0. This proves the first statement 
of Theorem 3.6. 

To prove the second statement, it is of course enough to show that for fixed 
p(a) > and all p(b) > sufficiently close to p(a), t{B 00 ) is bigger than one. 
By [83, Theorem 2.1-iv], we know that dimn(i?oo) can be forced arbitrarily close 
to one by taking p(h) sufficiently close to p{a). This, together with the following 
relation between Hausdorff dimension and thickness (see [59] for details) gives 
r(Soo) > 1: 

log 2 
-, r- < dim H (»). 

Finally, (9) follows since Boo is symmetric with respect to the origin in R (see 
proof of Theorem 3.4 above, and notice the appearance of E 2 in the definition of 

7). 



Let K^^n = Jui.n/ksT and h um = H Utn /kBT- Then the partition function ZL 



4.4 Proof of Theorem 3.8 

can be written as 

Z (N) = Tr (V W) f where T (N) = -Q Twj . (2Q) 

i=N 
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with 



% 



2 



exp ( Ku, 

exp (-K u>i ) 



cxp {-K Wji ) 
exp !«„,, 2 ' 



(21) 



(recall - we are using circular boundary conditions, that is, iST Wi jv+i = ^w,i and 
^w,jv+i = h u ,i). For a detailed derivation see, for example, [6]. Note that, un- 
like in the periodic case, the matrices {7L,i} do not commute and so cannot be 
diagonalized simultaneously. This is precisely what complicates the analysis, as 
the partition function, and hence the free energy function, cannot be computed 
explicitly. 

For a matrix M = iflij), define the one- norm and the infinity-norm, respec- 
tively, by 

ll M lli -=^2\ a ij\ and ll M lloo : =max{|o«|}- 



From (21) it is evident that for each 7L,i, HTL.iHoo = [7w,»]ii. Hence from (20) we 
obtain 



7-W 



From (22) and (12) we obtain 

1 



Nk B r 
Hence by Theorem 2.1 



log 



rW 



< Z,W < 



> F,W > 



rt") 



i 



Nk B r 



log 



7-W 



(22) 



(23) 



lim F,W 



JV- 



cxists as a function of temperature r € (0, oo) and is independent of u. From 
(22) and (23) it follows that F is strictly negative. It remains to prove that F is 
analytic in r. 

Since the limit in (23) is independent of the choice u €. Cl, we may choose 
uj s £ Q, where uj s is the special sequence as constructed in Section 2.1. From now 
on we drop to in all formulas, assuming implicitly that to = lu s . 

Since the limit in (23) exists as N — > oo, we may pass to a subsequence Nk = 
Fk, where Fk is the kth Fibonacci number, and by (20) it is enough to prove 
analyticity of the limit 



logTr(T (Ffc) ) . logTr(T (Ffc 
lim ; = lim )-, 



k— >oo 



F k 



k— >oo 



where <f> is the golden mean (the last equality follows since limfc_ i . 0O Fk/(f> 
Let us define 



(24) 
1). 



77 := -rTl, where di := \/det7i, and set 
di 



Xk := 



\i=F h I 
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Since % is unimodular, the trace map can now be applied to these matrices. For 
k > 3, we have (see, for example, [2] and references therein) 



Tr(T^)=(2f[ck\nof-Z( 



x 3l x 2 ,xi), (25) 



where -k is projection onto the first coordinate, and / is the Fibonacci trace map 
as defined in (14). 

Observe that the sequence {di} ieN is two- valued modulated by the sequence u, 
as defined in Section 2.1. We have 

log (2 ft*!*) log(2df- 1 ^- 2 ) j 1 
lira — = lim — = — log d\ -\ log di , 

k— ¥00 Fk k— yoo Ffc (p 1 + 

which is obviously analytic as a function of r. Thus in order to complete the proof, 
we need to prove that 

lOR 7T o f \X 3 Xo X~\ ) 

lim — - — - exists and is analytic in r. (26) 

k—^oo r ' fc 

The initial condition (x3,X2,Xi) defines an analytic curve (0,oo) -> I 3 as a 
map of temperature, r. Since the limit in (26) exists and is nonzero (by (23) 
and application of Theorem 2.1), {ir o / fc ~ 3 (x3,a;2,a;i)} fc>3 is unbounded, hence 
(x3,X2,xi) escapes under the action of /. 

Fix To € (0,oo). Consider the complexification of (x3,X2,Xi). Since for all 
ret, (xs,X2,Xi) escapes under the action of /, by Proposition 4.3, there exists 
a small neighborhood U around To in C, such that all z £ U have unbounded 
forward semi-orbit. Application of Proposition 4.4 gives (26). 

4.5 Proof of Theorem 3.10 

Notice that proof of the previous theorem gives a recursion relation for the partition 
function itself. Indeed, for k > 3, 



Z(*»= Uf[dAirof k - 3 (x 3 ,X2,x 1 ). 



Since we are interested in the zeros of Z, it is enough to consider 

2,(Fh) :=7r of k ~ 3 (x 3 ,X2,x 1 ) 

as a function of the variable 77, as defined in (13). In the new variables a, (3 and 
77 we have 

t-i, v / a/3fj a/3rj a 2 /3 2 - 1 \ 

/ (X 3 ,X2,X 1 )- 



2^^!' 2V^^T' V(« 4 -l)(/3 4 -l)/ ' 
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where 77 = 77 + 1). We took / 1 to obtain a simpler expression. Since the zeros of 
Z( F *' all lie on S 1 , we restrict the allowed values for 77 to [—2, 2]. 
Set 

(H = ( a ^ a ^ « 2 /3 2 ~ 1 \ / 27 s 

Notice that 7 is an analytic line in R 3 . For k > 3, let 

Z fc = {^e[-2,2]:7ro/ fc - 2 ( 7 (^))=0}. 

We are interested to know whether {Z^} converges in Hausdorff metric to some 
compact (nonempty) Z C [—2,2], and if so, what is the topological, measure- 
theoretic and fractal-dimensional structure of Z. 
Define 

Z = he [-2,2] : 0^(7(77)) is bounded! . 

Observe that Z may contain finitely many isolated points; let us call this set 
of isolated points T. We now prove that Zk > Z \ T in Hausdorff metric, 

k— >oo 

where T C T (we believe that T = 0, but we are currently unable to prove this 
rigorously). 

For the sake of convenience, define 

Z\=Z\f. (28) 

We also prove that Z satisfies i-iii of Theorem 3.10, which will complete the proof. 
Let us first prove that Z satisfies i-iii of Theorem 3.10. 

The value of the Fricke-Vogt invariant along 7(77) is given by 

(aW+4)(° a -/? 2 ) 2 ^ 

1{l{?i)) = 4(o4-l)(/34-l) • (29) 

Observe that 1(7(77)) = if and only if a = /3, in which case 7 lies entirely 
on Sq\ This corresponds to the case p(a) = p(b), i.e. the pure case in which the 
partition function and corresponding zeros can be computed explicitly. Otherwise, 
7(7(77)) > 0, i.e. 7 lies entirely in Uyxi'-'v- Moreover, 1(7(77)) is nonconstant as 
a function of 77. We are almost in the position where we could apply Theorem 4.5. 
Next, let us show that 7 cannot lie entirely on a center-stable manifold. Indeed, 
let us show that the endpoints of 7, that is, 7(— 2) and 7(2), escape to infinity. 
Indeed, at ±2 we have 

( ±a/3 ±a/3 o?f?-\ \ 

{x,y,z) :=7(±2) 



Clearly either \x\ > 1 or \y\ > 1. Using this, a direct computation shows that either 
fix, j/, z) or / _1 (x, y, z) satisfies that the first two coordinates in absolute value are 



21 



greater than one, and their product is greater than the third coordinate. It follows 
that these points escape to infinity (see the discussion in Section 4.2 following 
equation (19)). This also proves ii of Theorem 3.10, since rj = ±2 corresponds to 
?? = ±1. 

Next, observe that 

7(0)= ( 0,0, ; a2 P 2 - l ^\ ( 30 ) 

V V(« 4 -l)(/3 4 -l)/ 

is a periodic point of period six (the exact form of this expression is not important; 
in general points of the form where two of the coordinates are zero are periodic). It 
follows that 7 intersects center-stable manifolds at points other than the cndpoints. 
Finally, let us show that the intersection of 7 with the center-stable manifold 
containing the point 7(0) is transversal. For convenience, let us call this center- 
stable manifold W. From (29) we have 



' ) T 

7( 7 (0)) > and — 



= 0. 

7(0) 



Hence 7(0) G S v with V > and 7'(0) G T l(0] S v . On the other hand, intersec- 
tions of center-stable manifolds with surfaces S v for any V > are transversal 
(see Proposition 4.2 above). Therefore we need to prove that 7'(0) is transversal 
in T 7 ( )Sy to the tangent space at 7(0) of the smooth curve WT\S V . On the other 
hand, the curve Wr\S v is precisely the stable manifold to the hyperbolic periodic 
point 7(0) G S v . Since this periodic point is of period six, passing to / 6 , we have 
/ 6 (7(0)) = 7(0)) - a fixed hyperbolic point under / 6 . Hence the tangent space at 
7(0) to WDSy coincides with the eigenspace corresponding to the smallest (in ab- 
solute value) eigenvalue of Df®, Q s — since / is easily seen to be area preserving, at 
a hyperbolic fixed point the differential has three distinct eigenvalues: {1, A, a/ A}, 
where |A| > 1 and a G {±1}- As can easily be checked, the sought eigenvalue and 
the corresponding eigenvector are 



-4a 2 V4a 4 +T+8a 4 + l and 1,- 



V4a 4 + 1 + 2a 2 - 1 
2a _ ' 



with a = , = = 1 > 0. On the other hand, as is evident from (27), 7' is 
VC^-i K^-i) v ;. ; 

spanned by (1, -\/(/3 4 — l)(a 4 — l) -1 , 0. Since the second entry in this vector is 
strictly positive, while the second entry in the eigenvector is strictly negative, and 
all other entries coincide, the two vectors are clearly transversal in T^^Sy, and 
we're done. 

Now Theorem 4.5 can be applied and i of Theorem 3.10 follows. 

Next let us prove that Z satisfies iii of Theorem 3.10. Indeed, by Theorem 4.5 
it is enough to prove that for all p(a) and p(b) sufficiently close (i.e. for all a and 
f3 sufficiently close), 7 in (27) intersects the center-stable manifolds transversally. 
This follows from the same arguments as in [84, Proof of Theorem 2.5]. 
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Finally, let us prove that Zk > Z in Hausdorff metric. In place of /, let us 

k— >oo 

consider g := f G (the sixths iterate of /). That is, let us show that Z§k > Z 

k— ^oo 

in Hausdorff metric. This is done merely for the sake of technical convenience, as 
is apparent from the proof below, and does not result in loss of generality. 

Fix arbitrary e > and let U\, . . . , U mc be an open cover of Z with diam([/i) < 
e. It is enough to show that there exists A e G N such that for all k > K e , Z$k is 
contained in U£i Ui, and for every i G {1, . . . , m e }, Zek H Ui 7^ if Ui fl Z is not 
a finite set (i.e., if Ui does not only contain isolated points of Z). 

Every point of 7 n (U£i Ui) escapes to infinity. Hence there exists K\ G N 
such that for all k > K\, Z 6 k C U£i ^»- 

Consider the four singularities 

Pi = (1,1,1), P 2 = (-1,-1,1), P 3 = (1,-1,-1), P 4 = (-1,1,-1) 

of the Cay ley cubic Sq. The first is fixed under /, and the other three form a 
three-cycle. When the parameter V is turned on, the first singularity bifurcates 
into a hyperbolic periodic point of period two. The three cycle bifurcates into a 
hyperbolic six cycle. More precisely, there are four smooth curves in IR 3 that lie 
entirely in \J V>Q S V , one through each Pj, i = 1, ... ,4. Label these curves pi. 
Each pi is divided into two smooth branches by Pj; label these branches p\ and 
p\. Under /, p\ h-> p^, while the other six branches are cycled through a six-cycle. 
Points on these curves are hyperbolic. In particular, under g (recall: g = / ) these 
curves are fixed normally hyperbolic submanifolds of R 3 (for a brief discussion of 
normal hypcrbolicity, see Appendix B, and references therein for details). Each of 
p? r , i = 1, . . . ,4, intersects each Sy, V > 0, transversally and in only one point. 
This point of intersection is a fixed hyperbolic point for g (for more details see 
Appendix A.l). 

The stable manifolds to each of the six branches are members of the family 
of center-stable manifolds and are dense in this family (for details, see Appendix 
A.l below, and references therein). Moreover, among the six branches, there 
are those that lie below the plane {z = 0}, and those that lie above (for details, 
see the discussion on symmetries of the Fibonacci trace map in Appendix A. 2). 
Consequently, for each %q g {1, . . . , m e }, if Ui does not contain only isolated points 
of Z, then there is a pair of points p,q G Z n Ui al such that p lies on the stable 
manifold to a branch that lies below {z = 0}, and q lies on the stable manifold to 
a branch that lies above {z = 0}. Hence for all k sufficiently large, g k (p) lies below 
{z = 0}, and g k {q) lies above {z = 0}. Therefore, if a denotes the arc along 7 
that connects p and q (and a priori lies in Ui ), then for all such k, g(a) intersects 
{z = 0}. Hence for all such k, Z 6 k fl Ui ^ 0. Consequently, there exists K 2 G N 
such that for all k > K 2 , Z^ n Ui ^ for every i G {1, . . . , m e } such that Ui does 
not contain only isolated points of Z. 

Let K c = K\ + A'2 • This completes the proof. 
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Appendices 

A Normal hyperbolicity of six-cycles through sin- 
gularities of Sq and symmetries of / 

A.l Normal hyperbolicity of six-cycles through singularities 

of So 

As has been mentioned above, 5*0 contains four conic singularities; explicitly, they 
are 

Pi = (1,1,1), P 2 = (-1,-1,1), P 3 = (1,-1,-1), P 4 = (-1,1,-1). (31) 
The point Pi is fixed under T, while P 2 , P3 and P4 form a three cycle: 

Pl^+Pu P 2 ^ P 3 ^ P 4 ^ P2 

(which can be verified via direct computation). 

For each i G {i = 1, . . . ,4}, there is a smooth curve pi, containing no self- 
intersections, passing through the singularity p, such that pi \ Pi is a disjoint 
union of two smooth curves — call them p\ and p\ — with the following properties: 

• f(Pi) = Pi an d f{Pi) = fiPi)- I n particular, points of p{ r are periodic of 
period two, and p\ is fixed under /; 

• The six curves p/ r , i = 2, 3, 4, form a six cycle under /. In particular, points 

of Pi T , i = 2,1 
fixed under / 6 . 

The curve p\ is given explicitly by 



of pf , % = 2,3,4, are periodic of period six, and hence for i = 2,3,4, pi is 



Pi = \(x, ^—j, x) :xe (-00, 1/2) U (1/2, 00) J . (32) 

Expressions for the other three curves can be obtained from (32) using symmetries 
of / to be discussed below. 

It follows via simple computation that for any i = 1, . . . , 4 and any point p G 
p? T , the eigenvalue spectrum of Df£ is {1, A(p), 1/X(p)} with < |A(p)| < 1, where 
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Dfp denotes the differential of f 6 at the point p. The eigenspace corresponding to 
the eigenvalue 1 is tangent to pi at p. At Pj, the eigenvalue spectrum of Dfp is of 
the same form, and as above, the eigenspace corresponding to the unit eigenvalue is 
tangent to pi at Pj. It follows that the curves {pi} i=1 4 are normally hyperbolic 
one-dimensional submanifolds of R 3 , as defined in Section B.3. 

For V > 0, Sy f](pi U • • • U pi) consists of eight hyperbolic periodic points for 
/ which are fixed hyperbolic points for f 6 . The stable manifolds to these points 
form a dense sublamination of the stable lamination 

w s (n v )= |J w s ( x ), 

the union of global stable manifolds to points in Sly, the nonwandering set for / 
on Sy from Theorem 4.1. Moreover, this sublamination forms the boundary of the 
stable lamination. For details, see [14,15,20]. This extends to the center-stable 
lamination: the stable manifolds to the normally hyperbolic curves p\ , . . . , p± form 
a dense sublamination of the (two-dimensional) center-stable lamination and forms 
the boundary of this lamination. In particular, if 7 is a smooth curve intersecting a 
center-stable manifold, and this intersection is not isolated, then in an arbitrarily 
small neighborhood of this intersection, 7 intersects the stable manifolds of all 
eight curves: < pf > 

I J i— 1,...,4 

A. 2 Symmetries of / 

The following discussion is taken from our forthcoming paper [22]. 

Let us denote the group of symmetries of / 6 by S ymj an d the group of reversing 
symmetries of f e by Q rcv ; that is, 

G sym = {s e Diff(R 3 ) : s o / 6 o s- 1 = / 6 } , (33) 

and 

Srcv = {s G Diff (R 3 ) :sofo s- 1 = r 6 } , (34) 

where Diff (R 3 ) denotes the set of diffcomorphisms on R 3 . 
Observe that Q ICV ^ 0. Indeed, 

s(x,y,z) = (z,y,x) (35) 

is a reversing symmetry of /, and hence also of / 6 . Hence / 6 is smoothly conjugate 
to /~ 6 . It follows (see Appendix B) that forward-time dynamical properties of / 6 , 
as well as the geometry of dynamical invariants (such as stable manifolds) arc 
mapped smoothly and rigidly to those of /~ 6 . That is, forward-time dynamics of 
/ 6 is essentially the same as its backward-time dynamics. 

The group Q sym is also nonempty, and more importantly, it contains the fol- 
lowing diffcomorphisms: 

s 2 : (x,y,z) ^ (-x,-y,z), 

s 3 : (x, y, z) ^ (x, -y, -z), (36) 

s 4 : (x,y,z) i-» (-x,y,-z). 
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Notice that Si are rigid transformations. Also notice that 

s l {P 1 )=P l , 
and since Si is in fact a smooth conjugacy, we must have 

Si(pi) = pi- (37) 

For a more general and extensive discussion of symmetries and reversing symme- 
tries of trace maps, see [4]. 

Now from (37), (32) and conclusions of Section A.l it follows that for any V > 
0, there exist two points p,g€ ily fixed under / 6 , such that p lies below {z = 0}, 
q lies above {z = 0}, and W s (p) and W s {q) are dense in the stable lamination 

w s (n v ). 

B Background on uniform, partial and normal 
hyperbolicity 

B.l Properties of locally maximal hyperbolic sets 

A more detailed discussion can be found in [30,31,34-36]. 

A closed invariant set A C M of a diffcomorphism / : M — >• M of a smooth 
manifold M is called hyperbolic if for each x £ A, there exists the splitting T X A = 
Ex ® -^x invariant under the differential Df, and Df exponentially contracts 
vectors in E x and exponentially expands vectors in E%. The set A is called locally 
maximal if there exists a neighborhood U of A such that 

A=f]f l (U). (38) 

The set A is called transitive if it contains a dense orbit. It isn't hard to prove 
that the splitting E s x © E x depends continuously on x £ A, hence dim.(E x ' u ) is 
locally constant. If A is transitive, then dim(E x ,u ) is constant on A. We call the 
splitting E s x © Kjf a (k x ,k x ) splitting if dim(^ ,u ) = k s ' u , respectively. In case A 
is transitive, we'll simply write (fc s , k u ). 

Definition B.l. We call A C M a basic set for / e DifT(M), r > 1, if A is a 
locally maximal invariant transitive hyperbolic set for /. 

Suppose A is a basic set for / with (1, 1) splitting. Then the following holds. 

B.l.l Stability 

Let U be as in (38). Then there exists U C Diff (M) open, containing /, such 
that for all g £U, 

A, = fl 9 n (U) (39) 

n£Z 
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is (/-invariant transitive hyperbolic set; moreover, there exists a (unique) homeo- 
morphism H g : A — > A g such that 

H g of\ A =g\ Ag oH g . (40) 

Also H g can be taken arbitrarily close to the identity by taking U sufficiently small. 
In this case g is said to be conjugate to /, and H g is said to be the conjugacy. 

B.1.2 Stable and unstable invariant manifolds 

Let e > be small. For each x € A define the local stable and local unstable 
manifolds at x: 

W e s (x) = {ye M : d(f n (x), f"(y)) < e for all n > 0} , 

W?(x) ={y€M: d(f n (x), f n (y)) < e for all n < 0} . 

We sometimes do not specify e and write 

Wf oc (x) and W£ e (x) 

for W*{x) and W^(x), respectively, for (unspecified) small enough e > 0. For all 
ieA, W°£{x) is an embedded C r disc with T x W££{x) = E^ u . The global stable 
and global unstable manifolds 

W s (x) = |J r n (Wf oc (x)) and W u (x) = (J P{W? oc (x)) (41) 

nSN nGN 

are injectively immersed C r submanifolds of M. Define also the stable and unsta- 
ble sets of A: 

W S (A) = (J W s (x) and W U (A) = |J W(ic). (42) 

If A is compact, there exists e > such that for any x, y £ A, W/(a:) n ^"(y) 
consists of at most one point, and there exists 5 > such that whenever d(x, y) < S, 
x, y <G A, then W^{x) f~l W^y) 7^ 0. If in addition A is locally maximal, then 
W°{x)nW?(y) € A. 

The stable and unstable manifolds W^(x) depend continuously on x in the 
sense that there exists § s ' u : A — > Emb r (R, M) continuous, with $ s '"(a;) a neigh- 
borhood of x along W[q"(2;), where Emb r (R, M) is the set of C r embeddings of R 
into M [35, Theorem 3.2]. 

The manifolds also depend continuously on the dificomorphism in the following 
sense. For all g £ DifT(M) C r close to /, define $*'" : A g -> Emb r (R,Af) as we 
defined $ s,u above. Then define 

$£' u : A^Emb r (R,A/) 
by 

Then § s ' u depends continuously on 5 [35, Theorem 7.4]. 

27 



B.1.3 Fundamental domains 

Along every stable and unstable manifold, one can construct the so-called funda- 
mental domains as follows. Let W s (x) be the stable manifold at x. Let y £ W s {x). 
We call the arc 7 along W s (x) with endpoints y and f~ 1 (y) a fundamental domain. 
The following holds. 

• f(l) H W s (x) = y and f~ l {i) n W^(a-) = / -1 (y), and for any & G Z, if 
k < -1, then / fc ( 7 ) n W s (x) = 0; if fc > 1 then / fc (7) H W^(x) = iff x ^ y; 

• For any z £ W s (x), if for some k € N, / fe (z) lies on the arc along VF s (x) 
that connects 2 and y, then there exists n <G N, n < fc, such that f n {z) € 7. 

Similar results hold for the unstable manifolds. 

B.1.4 Invariant foliations 

A stable foliation for A is a foliation F s of a neighborhood of A such that 

1. for each x <G A, T(x), the leaf containing x, is tangent to E%.; 

2. for each x sufficiently close to A, f(J rs (x)) C J 7S (f(x)). 

An unstable foliation T u is defined similarly. 

For a locally maximal hyperbolic set AcM for / <G Diff (M), dim(M) = 2, 
stable and unstable C° foliations with C 1 leaves can be constructed; in case / € 
Diff (M), C 1 invariant foliations exist (see [59, Section A.l] and the references 
therein) . 

B.1.5 Local Hausdorff and box-counting dimensions 

For x € A and e > 0, consider the set W^ u n A. Its Hausdorff dimension is 
independent of x e A and e > 0. 
Let 

h s - u (A) = dim;/ (^'"(a;) n A). (43) 

For properly chosen e > 0, the sets W^' w (x) n A are dynamically defined Cantor 
sets, so 

h s - u (A) < 1 

(see [59, Chapter 4]). Moreover, h s ' u depends continuously on the diffcomorphism 
in the C 1 -topology [56]. In fact, when dim(M) = 2, these arc C r_1 functions of 
/ e Diff'(Af), for r > 2 [55]. 

Denote the box-counting dimension of a set T by dimBox(r). Then 

dim H (W°' u (x) n A) = dim B ox(W-7'"(20 n A) 

(see [56,73]). 
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B.2 Partial hyperbolicity 

For a more detailed discussion, see [32,60]. 

An invariant set A C M of a diffeomorphism / <E Diff r (M), r > 1, is called 
partially hyperbolic (in the narrow sense) if for each x € A there exists a splitting 
T X M = E^®E^,®E^ invariant under Df, and Df exponentially contracts vectors 
in E xl exponentially expands vectors in E x , and Df may contract or expand vec- 
tors from E x , but not as fast as in E* ,u . We call the splitting (k x , k x , k x ) splitting 
if dim(£'|' c ' u ) = k x ' c ' u , respectively. We'll write (k s ,k c ,k u ) if the dimension of 
subspaces does not depend on the point. 

B.3 Normal hyperbolicity 

For a more detailed discussion and proofs see [36] and also [60]. 

Let M be a smooth Ricmannian manifold, compact, connected and without 
boundary. Let / £ Diff r (M), r > 1. Let TV be a compact smooth suhmanifold 
of M, invariant under /. We call / normally hyperbolic on TV if / is partially 
hyperbolic on TV. That is, for each x E TV, 

T X M = E s x @E c x ®E x l 

with E x = T X TV. Here E x ,c,u is as in Section B.2. Hence for each x € TV one 
can construct local stable and unstable manifolds W"{x) and W™(x), respectively, 
such that 

1. seW£ c (s)nWEc(a;); 

2. T x Wf oc (x) = E°(x), T a W$Jp) = £»(*); 

3. for n > 0, 

d(f n (x),f n (y)) > for all y € Wf oc {x), 

n—too 

d(f- n (x), /"»(!/)) ► for all j, e W? oc (x). 

n—^oo 

(For the proof see [60, Theorem 4.3]). These can then be extended globally by 

W'(x) = |J f- n f(W{ oc (x)); 

nGN 

W u (x) = |J r(W£ c (x)). 



The manifold VF s (x) is referred to as the strong-stable manifold, while VF u (x) 
is called the strong-unstable manifold; sometimes to emphasize the point x, we 
add at x. 

Set 

CcW = (J W£c(s) and ^locW = (J WiocO*). 

xGJV KEN 
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Theorem B.2 (Hirsch, Pugh and Shub [36]). The sets W{£.(N) and W{%.(N), 

restricted to a neighborhood of N , are smooth submanifolds of M . Moreover, 

1. WjQ S c (iV) is f -invariant and Wj™ is f~ x -invariant; 

2. N = W 1 Z(N)f)W&(N); 

3. For every x&N, T x W^ c cn {N) = E s x > u © T X N ; 

4- Wfo c (N) (W{^.(N)) is the only f -invariant (f^ 1 -invariant) set in a neigh- 
borhood of N ; 

5. W™ C (N) (respectively, W™ C (N)) consists precisely of those points y G M 
such that for all n > (respectively, n < 0), d(f n (x),f n (y)) < e for some 
e>0. 

6. W^ CU (N) is foliated by {W^(x)} x&N . 
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